#ifndef RASTERTOFIELD_HPP
#define RASTERTOFIELD_HPP

#include "meshTypes_fwd.hpp"
#include "raster.hpp"

namespace hms::raster {

scalar rasterValueAtCoords(
	const Raster& raster,
	const Ref<const Array2s>& coords
);


template<typename MeshType>
ArrayXXs rasterToFieldValues(
	const Raster& raster,
	const MeshBase<MeshType>& mesh
);

} // namespace hms::raster

#endif